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Abstract 

This study presents the distribution of bacterial and fungal abundances in long-term copper (Cu) con- 
taminated soils in Topolnitsa-Pirdop valley — a highly industrialized zone with a number of mines and 
processing plants for copper and other non-ferrous metals. The bacterial (16S rRNA gene copies) and fun- 
gal (ITS rRNA gene copies) were estimated using quantitative PCR technique in five topsoils, differently 
Cu contaminated (ranging from 28.05 to 198.9 mg kg"). Bacterial abundance varied in a range of 1.68 
x 10" to 3.24 x 10!'16S rRNA genes, whereas fungi amounted from 1.95 x 10% to 6.71 x 10° ITS rRNA 
genes. Fungal and bacterial abundances were significantly (fungi) and insignificantly (bacteria) influenced 
by Cu contamination. The fungal/bacterial ratio related negatively with soil Cu, which shifted microbial 
communities’ structure towards bacterial dominance. Since the ratio between bacteria and fungi are vital 
in explaining many soil functions, the calculated changes in this ratio indicated deterioration in soil qual- 


ity, being of primary importance for plant production. 
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Introduction 


Heavy metals (HMs) are natural constituents of the environment, but intensive use 
for human purposes has altered their geochemical cycles and biochemical balance. 
The extraction and processing of mineral raw materials are one of the main sources 
of soil and water pollution with HMs and other toxic elements. In this respect, of 
particular interest are the open-pit mines and their generated waste materials, which 
cause serious environmental damage and significant changes to the landscapes of 
the impacted regions (Dabeva et al. 2012). Prolonged exposures to HMs in high 
concentrations are harmful to microorganisms, plants and animals. Furthermore, 
HMs can be absorbed by food crops, and entering into the food chain has proven 
to be a potential health hazard for plants and humans (Huang et al. 2018; Ali 
et al. 2019; Afonne and Ifediba 2020; Hasan et al. 2020). Microorganisms are 
highly diverse and ubiquitous in soil ecosystems and participate in a variety of 
key ecosystem functions such as nutrient cycling, structuring of soil aggregates 
and biomass production. A stable microbial community contributes essentially to 
stabilizing soil structure and maintaining soil ecosystem services (Bissett et al. 2013). 
Studies have shown that microorganisms can significantly promote the circulation 
of soil nutrients, maintain soil fertility, and improve crop health (Fierer 2017). The 
negative impacts of HMs on soil microbial communities have been reported in 
several studies. HMs have been shown to be harmful to soil microorganisms and 
may lead to significant changes in microbial diversity (Chodak et al. 2013; Zampieri 
et al. 2016), richness (Cui et al. 2018), abundance (Deng et al. 2015), metabolic 
activity (Chen et al. 2014; Hong et al. 2015; Zampieri et al. 2016; Fajardo et al. 
2019) and structure of the communities (Cui et al. 2018; Feng et al. 2018; Zhang 
et al. 2018). 

Copper is one of the most common soil contaminants (Griffiths and Philippot 
2013) among heavy metals. Soils represent the largest sink of Cu, being released to the 
environment by anthropogenic activities, such as sewage irrigation, mining activities, 
municipal waste disposal, and intensive use of pesticides and herbicides (Smith 2009; 
Zhuang et al. 2009; Ballabio et al. 2018). A moderate level of Cu is regarded as a mi- 
cronutrient and it is essential for microorganisms to carry out their metabolic activities 
(Giller et al. 2009). However, excessive inputs of Cu to soil ecosystems could persist 
for a long time after their introduction and cause toxic effects on soil microorganisms 
(Giller et al. 2009). Some authors found that Cu contamination had significant effects 
on bacterial abundance and diversity (Li et al. 2015), as well as fungal diversity (Zhang 
et al. 2022). Therefore, understanding the responses of soil microbial communities to 
Cu contamination is essential to counteract its negative effect on ecosystem functions 
and services. 

The aim of this study was to assess the distribution of bacterial and fungal abun- 
dances in Cu long-term contaminated soils from the valley of Topolnitsa-Pirdop and 
its relationship with the local soil properties. 
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Materials and methods 


Study area and soil sampling 


Topolnitsa River and its tributaries (Medetska, Zlatishka, Pirdopska and Bunovska 
Rivers) were influenced by the ore industry, agriculture activities, and domestic waste- 
water discharge, which led to the deterioration of its ecological status according to the 
River Basin Management Plan (2016-2021) (East Aegean River Basin Directorate. 
MoEW). Three large metal mining and metal processing enterprises are located in this 
region- “Elatsite Med”, “Dundee Precious Metals” and “Aurubis”. Additionally, in the 
area, the open-pit “Medet” is located, and although the mine was closed in 1992, it is 
considered the main source of pollution in this area. 

The study area is located in the Topolnitsa-Pirdop valley (Fig. 1). Topsoil samples 
(0-20 cm) with different Cu concentrations (ranging from 28.05 to 198.9 mg ke") 
were collected in May 2021. Five sampling sites have been selected from a monitor- 
ing map (Kancheva et al. 2018), considering the sources of Cu pollution as follows: 
S1 — located in the valley of Medetska River at 710 m from the “Medet” mine tailing 
(42°39'12.9"N, 24°09'04.2"E), S2 —located in the valley of Topolnitsa River (after Zlat- 
ishka River) at 250 m from the “Medet” mine tailing (42°39'31.8"N, 24°08'26.6"E), 
S3 — located in the valley of Pirdopska River at 750 m from the “Medet” mine tailing 
(42°39'48.2"N, 24°08'11.3E), S4—located in the valley of Bunovska River at 2200 m 
from the mine tailings of “Elatsite-Med” A.D. (42°36'47.5"N, 24°00'44.8"E), and 
S5-located in the valley of Topolnitsa River-next to Muhovo village at 25-30 km 
straight from S4 and S1 (42°24'45.5"N, 23°59'54.4"E). Samples were taken in three 
replicates per site. All sampled soils were classified as Fluvisols (FAO 2015). 


Soil physicochemical properties and copper concentration 


Soil pH was potentiometrically measured in 0.1 M CaCl, soil suspension according to 
ISO 10390:2005. The soil organic carbon (SOC) was determined by sulfochromic oxi- 
dation according to BSS ISO 14235:2002. Nitrate (NO,-N) and ammonium (NH,-N) 
nitrogen were determined according to Keeney and Nelson (1982). Inorganic phosphate 
(P,O,) was determined by the method of Olsen (1982). Water content (WtC) was calcu- 
lated after oven drying (105 °C). The concentration of Cu was measured by inductively 
coupled plasma-atomic emission spectrometry (ICP-AES) according to ISO 22036:2008. 


DNA extraction 


The genomic DNA was extracted from 0.5 g soil, using the E.Z.N.A soil DNA kit 
(Omega Bio-tek, USA) according to the manufacturer‘s recommended protocol. The 
quality and quantity of the extracted DNA were subsequently assessed by Qubit 4 
fluorometer (Invitrogen) and 1% agarose gel electrophoresis. 
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gend 

$1 - Medetska River (outfall) 

$2 - Topolnitsa River - after Zlatishka River outfall 
S3 - Pirdopska River (outfall) 

$4 - Bunovska River (outfall) 

$5 - Topolnitsa River - next to Muhovo village 


Figure |. Map of the study area and sampling sites. 


Real-time Q-PCR 


Q-PCR analyses were performed using iTaql@M Universal SYBRGreen Supermix 
(BioRad) on Rotor Gene 6000 (Corbett Life Science) according to the manufactur- 
er’s protocols. Bacterial abundance was estimated using 0.4 uM of primers Eub338f 
(5’-ACTCCTACGGGAGGCAGCAG-3’) and Eub518r (5’-ATTACCGCGGCT- 
GCTGG-3’) (Aleksova et al. 2020). The three-step thermal program was used at 
95 °C for 5 min,45 cycles of denaturing at 95 °C for 10 s, annealing at 61 °C for 30 
s, extension at 72 °C for 30 s, and final extension at 72 °C for 10 min. The plasmid 
DNA from the bacterial standard of Uncultured « — proteobacterium clone BuhD-104 
(FM877535.1) was used to generate a linear calibration curve of the threshold cycle 
versus the gene copy numbers of four-point serial decimal dilution — 10°, 10°, 10° 
and 10%. All measurements were run in triplicates. Data were expressed as gene copy 
numbers per gram of dry soil, and the specificity of amplicons was confirmed by melt- 
ing curve and agarose gel electrophoresis. Amplification efciency was estimated to be 
109% (R’°=0.99). The results were processed via ROTOR-GENE Q SERIES, Software 
version 2.3.1. 

The fungal abundance was estimated using 0.4 uM primers ITS1 (5’°-TCCGTAG- 
GTGAACCTGCGG-3’) and ITS4 (5’°-TCCTCCGCCTTATTGATATGC-3’) (White 
et al. 1990). The following thermal program was used: 95 °C for 5 min, 45 cycles of 95 
°C for 30 s, annealing 55 °C for 30 s, 72 °C for 60 s, and final extension at 72 °C for 10 


min. 
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A four-point serial decimal dilution ofa plasmid DNA of Uncultured Basidiomycota 
clone LS_Az0O_D1_31 (M1T785782.1) was used as a standard to generate a linear cali- 
bration curve of the threshold cycle versus a number of gene copies ranging from 10° to 
10°. All measurements were run in triplicates. Data were expressed as gene copy num- 
bers per gram of dry soil. Amplification efficiency was estimated to be 92% (R’*=0.99). 
The results were analysed via ROTOR-GENE Q SERIES, Software version 2.3.1. 

The relative ratio of fungi/bacteria was calculated as the ratio of the fungal gene 
copy numbers to bacterial gene copy numbers (Fierer et al. 2005). 


Statistical analysis 


Pearson correlation analysis was used to find the relationships between soil metrics. Hierarchi- 
cal clustering (Algorithm: UPGMA; Similarity index: Bray-Curtis) was conducted to evalu- 
ate the distance between soil physical environments. Non-metric multidimensional scale 
(nMDS) ordination was applied to assess the similarity between abundance-weighted mi- 
crobial communities. Statistical analyses were performed using the software PAST (p<0.05). 


Results 


Soil environmental variables 


The values of studied soil properties are shown in Table 1. The soil pH was moderately 
acidic — pH 5.57 (S1 and S2) to 5.76 (S4). SOC ranged from 1.22 g kg" (S3) to 2.67 g 
ke! (S2). Concentrations of inorganic nitrate (except $3) and ammonium (except S2 and 
S4) ions and inorganic phosphates were relatively equally distributed among the soils. 
Cu concentrations varied among the soils as follows: the concentrations in $1 and 
S3 exceeded the maximum permissible concentration (MPC) of 80 mg kg" under Bul- 
garian Regulation 3/2008 (2008) (http://eea.governement.bg/bg/legislation/soil). The 
most contaminated soils are located near the tailing of mine “Medet” in the valleys of 
the Medetska and Pirdopska Rivers. In the other soils, S2, S5 and S4, the concentrations 
of Cu were under or around MPC. In this regard, $5 can be assumed as a reference site. 


Table 1. Soil physicochemical properties and Cu concentrations in the region of Topolnitsa River and 


its tributaries. 


Soil sample Soil physicochemical properties 

pH WtC (%) SOC gkg' NO,-N mg kg’ NH,-N mg kg! P,O, mgkg' Cu mg kg" 
$1 D207 13.00 2.56 11.44 7.82 ie Fe) 198.90 
S3 5.60 11.67 12. 0.65 8.32 ae A 110.00 
S4 5.76 12.33 2.16 8.84 5.74 3.40 82.40 
52 5sid7 12233 2.67 11.12 3.07 3.47 67.40 
$5 5.61 10.00 2.05 tND 10.74 2.78 28.50 


+ND — No data. 
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Pearson correlation analysis indicated a linear relationship (-0.61; p=0.015) be- 
tween the level of soil Cu contamination and the distance to the source of pollution. 

A cluster analysis of soil abiotic factors was conducted to visualize the similarity 
among the tested soils (Fig. 2). 

Two main clusters were demonstrated: cluster I containing soil S5 with the lowest 
Cu concentration and cluster II, consisting of other soils sub-clustered according to the 
Cu threshold concentration of 100 mg kg’. This pattern of soil clustering identified 
the major role of Cu in structuring soil environments. 


Soil bacterial abundance 


Bacterial abundance varied in a relatively narrow range. The highest number was ob- 
served in S4 (3.24 x 10"'), and the lowest one in S5 (1.68 x 10"') (Fig. 3). Pearson 
correlation analysis did not indicate a significant correlation (0.31; p=0.26) between 
bacterial abundance and soil Cu concentrations. 


Soil fungal abundance 


The highest fungal abundance was detected in uncontaminated soil $5 (6.71 x 10°) 
and the lowest one in $3 (1.95 x 10%) (Fig. 4). The fungal gene copy number variance 
among the sampling soils was higher (Coef. variance — 40.48) than that of bacterial 
gene copy number (Coef. variance — 25.07), indicating possible resistance of bacteria 
to Cu toxicity. 
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Figure 2. Cluster analysis of soil physical environments. 
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Figure 3. A bacterial abundance of 16S rRNA gene. Error bars are the standard deviations of the mean 


for the three replicates. 
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Figure 4. An abundance of ITS rRNA gene. Error bars are the standard deviations of the mean for the 


three replicates. 


Pearson correlation analysis showed a significant negative correlation (-0.76; 
p=0.001) between fungal abundance and soil Cu contamination. 

The highest relative fungal to bacterial abundance was observed in the uncon- 
taminated soil S5 (3.98 x 10°), and the lowest ratio was recorded in the highly con- 
taminated soil S3 (6.80 x 10%) (Fig. 5). Fungal to bacterial ratio related negatively 
with soil Cu (-0.66; p=0.007), probably as a result of the higher Cu impact on the 


fungal abundance. 
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mean for the three replicates. 


Relationships between microbial abundance, soil properties and Cu 
concentrations 


Statistical non-metric multidimensional scaling (nMDS) was performed on patterns 
of fungal and bacterial abundances to assess soil similarities with respect to their mi- 
crobiology (Fig. 6). 

Plots separated the soils into four distinct abundance-weighted microbial commu- 
nities The soils S2 and S4 grouped closely, indicating similar microbial abundances. 
The other soils were placed individually both away from the other and from the group 
S2 and S4. This distribution was probably due to differences in microbial abundance- 
weighted communities. Our results showed that the important factors influencing soil 
microbial abundance were soil Cu concentration and water content in S1, and soil 
organic carbon and phosphates concentrations in $2 and S4 (Fig. 6). 


Discussion 


The present study was focused on the distribution of microbial abundance asa response to the 
influence of Cu contamination in the area of Topolnitsa-Pirdop valley — a highly industrial- 
ized zone with a number of mines and processing plants for Cu and other non-ferrous metals. 

The soil concentrations of the main contaminant Cu varied from 28.5 mg kg" to 
189.9 mg ke", being above MPC in three (S1, S3, and S4) of five tested soils. The mode 
of clustering highlighted the significant effects of Cu on soil environments, grouping 
soils in clusters or sub-clusters depending on the level of contamination. Additionally, 
soil pH (acidic) may influence soil microbial communities by primarily suppressing 
the bacterial diversity (Fierer and Jackson 2006) and increasing the mobility/toxicity 
(Wang et al. 2022) of HMs, including Cu. We expected also a high manifestation of 
Cu toxicity, taking into account the low levels of SOC (Chen et al. 2014). 
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Figure 6. Microbial abundances in soils of Topolnitsa-Pirdop valley as indicated by non-metric multi- 
dimensional scaling plots. 


The results showed high bacterial abundance in soils, ranging from 1.68 x 10'' to 
3.24 x 10'' 16S rRNA gene copy numbers, being one-fold higher compared to the 
abundances estimated in our previous studies on mining activities in the area, where the 
soil Cu concentrations ranged from 53 mg kg to 860 mg keg" (Aleksova et al. 2020; 
Palov et al. 2020; Nikolova et al. 2022). We suggested that this difference resulted from 
the ranges of soil Cu concentrations, but also microbial abundances might be edaphic- 
and climatic-dependent. In this respect, our results confirmed the findings of Yin et al. 
(2015), that long-term contamination with high Cu concentrations can decrease soil 
bacterial abundance. In our study, Pearson correlation analysis showed a lack of sig- 
nificant correlation between bacterial abundance and soil Cu concentration, probably 
due to the relatively low levels and long history of soil contamination. Hence, we could 
assume that bacteria developed mechanisms for adaptation to the newly created soil en- 
vironments, and relatively low Cu concentrations can even stimulate their proliferation 
(positive correlation between Cu and relative bacterial abundance) (Zhao et al. 2019). 

Similar to other authors (Chen et al. 2014), we found that the fungal abundance 
was sensitive to heavy metal contamination stress in the opposite to bacteria, proved 
by the significant and negative relationship between fungi and Cu. We suggested that 
the different sensitivity of bacteria and fungi to Cu resulted from their patterns of 
microbial- and HM- distributions among soil fractions. Several authors reported that 
bacteria dominated in silt and clay fractions, whereas fungi inhabited mainly the coarse 
sand fraction (Kandeler et al. 2000; Sessitsch et al. 2001; Poll et al. 2003). Using the 
Nemerow pollution index in their study, Chen et al. (2014) concluded that HMs (Cu, 
Pb, Zn) and their bioavailable forms were distributed in the order: coarse sand > clay 
> silt > fine sand. The similar mode of fungi and Cu distribution among soil fractions 
might be interpreted as a prerequisite for the higher Cu toxicity compared to bacteria. 
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The higher exerted Cu toxicity on fungi reflected also on the fungal/bacterial ratio, 
which decreased with increasing soil Cu contamination. According to some authors 
(Rajapaksha et al. 2004; Chen et al. 2014), the fungal/bacterial ratio is a sensitive indi- 
cator for soil health and its decrease in Cu-contaminated soils is a sign of deteriorated 
soil quality. Chen et al. (2014) reported a similar pattern of decreasing fungal/bacterial 
ratios under HM contamination in their study. 

Although microbial communities were differently abundance-weighted (except 
that of S2 and S4) according to the ordination procedure nMDS, it demonstrated that 
only the highest soil Cu concentration (189.9 mg kg’) influenced significantly soil 
microbial communities (S1). The abundances of soil microorganisms from the other 
sampling sites were under the influence of local soil abiotic factors. 


Conclusion 


The present study showed that microbial abundance, especially fungal, was signifi- 
cantly affected by long-term Cu contamination of Fluvisols. The Cu shifted microbial 
communities’ structure towards bacteria, suggesting that in this case, bacteria could be 
better at developing Cu resistance than fungi. Further studies should be implemented 
to clarify microbial functional responses to Cu. 


Acknowledgements 


This study was supported by the National Research Fund of the Bulgarian Ministry of 
Education and Science under Grant [KP-06-M41/3]. 


References 


Afonne OJ, Ifediba EC (2020) Heavy metals risks in plant foods — need to step up precau- 
tionary measures. Current Opinion in Toxicology 22: 1-6. https://doi.org/10.1016/j.co- 
tox.2019.12.006 

Aleksova M, Palov D, Dinev N, Boteva S, Kenarova A, Dimitrov R, Radeva G (2020) Bacterial 
abundance along a gradient of heavy metal contaminated soils in the region of Zlatitsa- 
Pirdop valley, Western Bulgaria. Dokladi na Bulgarskata Akademia na Naukite 73(3): 
433-442. https://doi.org/10.7546/CRABS.2020.03.18 

Ali H, Khan E, Ilahi I (2019) Environmental Chemistry and Ecotoxicology of Hazardous 
Heavy Metals. Environmental Persistence, Toxicity, and Bioaccumulation. Journal of 
Chemistry 2019(4): 1-14. https://doi.org/10.1155/2019/6730305 

Ballabio C, Panagos P, Lugato E, Huang JH, Orgiazzi A, Jones A, Fernandez-Ugalde O, Borrelli 
P, Montanarella L (2018) Copper distribution in European topsoils: An assessment based 
on LUCAS soil survey. The Science of the Total Environment 636: 282-298. https://doi. 
org/10.1016/j.scitotenv.2018.04.268 


Distribution of microbial abundance 33 


Bissett A, Brown MY, Siciliano SD, Thrall PH (2013) Microbial community responses to 
anthropogenically induced environmental change: Towards a systems approach. Ecology 
Letters 16(Suppl 1): 128-139. https://doi.org/10.1111/ele.12109 

BSS ISO 14235 (2002) Soil quality. Organic carbon determination by sulfo-chrome oxidation. 
https://bds-bg.org/en/project/show/bds:proj:30752 

Bulgarian Regulation 3/2008 (2008) On the norms for the permissible content of harmful 
substances in Bulgarian soils. https://eea.government.bg/ 

Chen Z, He EK Zhang X, Sun X, Zheng J, Zheng J (2014) Heavy metal pollution decreases 
microbial abundance, diversity and activity within particle-size fractions of a paddy soil. 
FEMS Microbiology Ecology 87(1): 164-181. https://doi.org/10.1111/1574-6941.12212 

Chodak M, Gotebiewski M, Morawska-Ploskonka J, Kuduk K, Niklinska M (2013) Diversity 
of microorganisms from forest soils differently polluted with heavy metals. Applied Soil 
Ecology 64: 7-14. https://doi.org/10.1016/j.apsoil.2012.11.004 

Cui H, Liu LL, Dai JR, Yu XN, Guo X, Yi SJ, Zhou DY, Guo WH, Du N (2018) Bacterial com- 
munity shaped by heavy metals and contributing to health risks in cornfields. Ecotoxicology 
and Environmental Safety 166: 259-269. https://doi.org/10.1016/j.ecoenv.2018.09.096 

Dabeva V, Georgieva D, Petrova S (2012) Investigation of the impact of contaminated water 
from the Topolnitsa River catchment on some agricultural crops. Ekologiyata — nachin na 
mislene 4: 9-10. [in Bulgarian] 

DenguL,~Zene G,-Fan: Geluel,” Chenkx,/ Ghent VM Wut Jie. XG: He ¥+@015) Response 
of rhizosphere microbial community structure and diversity to heavy metal co-pollution 
in arable soil. Applied Microbiology and Biotechnology 99(19): 8259-8269. https://doi. 
org/10.1007/s00253-015-6662-6 

Fajardo C, Costa G, Nande M, Botias P, Garcia-Cantalejo J, Martin M (2019) Pb, Cd, and 
Zn soil contamination: Monitoring functional and structural impacts on the microbiome. 
Applied Soil Ecology 135: 56-64. https://doi.org/10.1016/j.apsoil.2018.10.022 

FAO (2015) World reference base for soil resources 2014. International soil classification system 
for naming soils and creating legends for soil maps. Update 2015. World Soil Resources 
Reports No. 106. Food and Agriculture Organization of the United Nations, Rome, Italy, 
192 pp. http://www.fao.org/3/i3794en/13794en.pdf 

Feng M, Adams JM, Fan KK, Shi Y, Sun R, Chu H (2018) Long-term fertilization influ- 
ences community assembly processes of soil diazotrophs. Soil Biology & Biochemistry 126: 
151-158. https://doi.org/10.1016/j.soilbio.2018.08.021 

Fierer N (2017) Embracing the unknown: Disentangling the complexities of the soil micro- 
biome. Nature Reviews. Microbiology 15(10): 579-590. https://doi.org/10.1038/nrmi- 
ero2 O78 7 

Fierer N, Jackson RB (2006) ‘The diversity and biogeography of soil bacterial communities. 
Proceedings of the National Academy of Sciences of the United States of America 103(3): 
626-631. https://doi.org/10.1073/pnas.0507535103 

Fierer N, Jackson JA, Vilgalys R, Jackson RB (2005) Assessment of soil microbial community 
structure by use of taxon-specific quantitative PCR assays. Applied and Environmental 
Microbiology 71(7): 4117-4120. https://doi.org/10.1128/AEM.71.7.4117-4120.2005 

Giller KE, Witter E, Steve PR. McGrath SP (2009) Heavy metals and soil microbes. Soil Biology 
& Biochemistry 41(10): 2031-2037. https://doi.org/10.1016/j.soilbio.2009.04.026 


34 Michaella Petkova et al. / BioRisk 20: 23-35 (2023) 


Griffiths BS, Philippot L (2013) Insights into the resistance and resilience of the soil micro- 
bial community. FEMS Microbiology Reviews 37(2): 112-129. https://doi.org/10.1111/ 
j.1574-6976.2012.00343.x 

Hasan AB, Reza AHM, Kabir S, Siddique MAB, Ahsan MA, Akbor MA (2020) Accumula- 
tion and distribution of heavy metals in soil and food crops around the ship breaking area 
in southern Bangladesh and associated health risk assessment. SN Applied Sciences 2(2): 
1-18. https://doi.org/10.1007/s42452-019-1933-y 

Hong C, SiY, Xing Y, Yang L (2015) Illumina MiSeq sequencing investigation on the contrasting 
soil bacterial community structures in different iron mining areas. Environmental Science 
and Pollution Research International 22(14): 10788-10799. https://doi.org/10.1007/ 
s11356-015-4186-3 

Huang Y, Chen Q, Deng M, Japenga J, Li T, Yang X, He Z (2018) Heavy metal pollution 
and health risk assessment of agricultural soils in a typical peri-urban area in southeast 
China. Journal of Environmental Management 207: 159-168. https://doi.org/10.1016/j. 
jenvman.2017.10.072 

ISO 22036 (2008) Soil quality — Determination of trace elements in extracts f soil by induc- 
tively coupled plasma — atomic emission spectrometry (ICP-AES). https://www.iso.org/ 
standard/40653.html 

Kancheva V, Dinev N, Gadjalska N (2018) Risk assessment of copper contamination in 
soil and surface water in Topolnitsa river catchment. Bulgarian Journal of Soil Science. 
Argochemistry and Ecology 52(2): 3-9. [in Bulgarian] 

Kandeler E, Tscherko D, Bruce KD, Stemmer M, Hobbs PJ, Bardgett RD, Amelung W (2000) Struc- 
ture and function of the soil microbial community in microhabitats of a heavy metal polluted 
soil. Biology and Fertility of Soils 32(5): 390-400. https://doi.org/10.1007/s003740000268 

Keeney DR, Nelson DW (1982) Nitrogen — inorganic forms. In: Page AL (Ed.) Methods 
of Soil Analysis: Part 2 Chemical and Microbiological Properties. American Society of 
Agronomy, Soil Science Society of America, Madison, 643-698. https://doi.org/10.2134/ 
agronmonogr9.2.2ed.c33 

Li J, Ma Y, Hu H-W, Wang JT, Liu Y-R, He Ji-Z (2015) Field-based evidence for consistent re- 
sponses of bacterial communities to copper contamination in two contrasting agricultural 
soils. Frontiers Microbiology 6(31): 1-11. https://doi.org/10.3389/fmicb.2015.00031 

Nikolova R, Petkova M, Dinev N, Kenarova A, Boteva S, Berov D, Radeva G (2022) Correlation 
between bacterial abundance, soil properties and heavy metal contamination in the area 
of non-ferrous metal processing plant, Southern Bulgaria. BioRisk 17: 19-30. https://doi. 
org/10.3897/biorisk.17.77458 

Olsen SR (1982) Phosphorus. In: Page AL (Ed.) Methods of Soil Analysis: Part 2 Chemical 
and Microbiological Properties. American Society of Agronomy, Soil Science Society of 
America, Madison, 403-430. https://doi.org/10.2134/agronmonogr9.2.2ed.c24 

Palov DD, Aleksova MR, Nikolova RN, Dinev NS, Kenarova AE, Boteva SB, Radeva GS 
(2020) Relationships between Soil Microbial Activity, Bacterial Diversity and Abiotic 
Factors Along the Heavy Metal Contamination Gradient. Ecologia Balkanica 3: 31-39. 

Poll C, Thiede A, Wermbter N, Sessitsch A, Kandeler E (2003) Micro-scale distribution of mi- 


croorganisms and microbial enzyme activities in a soil with long-term organic amendment. 


Distribution of microbial abundance 33 


European Journal of Soil Science 54(4): 715-724. https://doi.org/10.1046/j.1351- 
0754.2003.0569.x 

Rajapaksha RMCP, Tobor-Kaplon MA, Baath E (2004) Metal Toxicity Affects Fungal and 
Bacterial Activities in Soil Differently. Applied and Environmental Microbiology 70(5): 
2966-2973. https://doi.org/10.1128/AEM.70.5.2966-2973.2004 

River Basin Management Plan (2016-2021) East Aegean River Basin Directorate. MoEW. https:// 
www.moew.government.bg/en/water/management-pans/river-basin-management-plans/ 

Sessitsch A, Weilharter A, Gerzabek MH, Kirchmann H, Kandeler E (2001) Microbial popu- 
lation structures in soil particle size fractions of a long-term fertilizer field experiment. 
Applied and Environmental Microbiology 67(9): 4215-4224. https://doi.org/10.1128/ 
AEM..67.9.4215-4224.2001 

Smith SR (2009) A critical review of the bioavailability and impacts of heavy metals in munici- 
pal solid waste composts compared to sewage sludge. Environment International 35(1): 
142-156. https://doi.org/10.1016/j.envint.2008.06.009 

Wang YF, Li X, Datta R, Chen J, Du Y, Du DL (2022) Key factors shaping prokaryotic com- 
munities in subtropical forest soils. Applied Soil Ecology 169: e104162. https://doi. 
org/10.1016/j.apsoil.2021.104162 

White T, Bruns T, Lee S, Taylor J (1990) Amplification and direct sequencing of fungal riboso- 
mal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ (Eds) 
PCR protocols: A Guide to Methods and Applications, Academic Press, Inc., New York, 
315-322. https://doi.org/10.1016/B978-0-12-372180-8.50042-1 

Yin H, Niu J, Ren Y, Cong J, Zhang X, Fan FE, Xiao Y, Zhang X, Deng J, Xie M, He Z, Zhou J, 
Liang Y, Liu X (2015) An integrated insight into the response of sedimentary microbial 
communities to heavy metal contamination. Scientific Reports 5(1): e14266. https://doi. 
org/10.1038/srep14266 

Zampieri BDB, Pinto AB, Schultz L, de Oliveira MA, de Oliveira AJ (2016) Diversity and 
Distribution of Heavy Metal-Resistant Bacteria in Polluted Sediments of the Araca Bay, Sao 
Sebastiao (SP), and the Relationship Between Heavy Metals and Organic Matter Concen- 
trations. Microbial Ecology 72(3): 582-594. https://doi.org/10.1007/s00248-016-0821-x 

Zhang X, Johnston ER, Barberan A, Ren Y, Wang Z, Han X (2018) Effect of intermediate distur- 
bance on soil microbial functional diversity depends on the amount of effective resources. En- 
vironmental Microbiology 20(10): 3862-3875. https://doi.org/10.1111/1462-2920.14407 

Zhang X, Fu G, Xing S, Fu W, Liu X, Wu H, Zhou X, Ma Y, Zhang X, Chen B (2022) Structure 
and diversity of fungal communities in long-term copper-contaminated agricultural soil. 
The Science of the Total Environment 806(3): 151302. https://doi.org/10.1016/j.scito- 
tenv.2021.151302 

Zhao X, Huang J, Lu J, Sun Y (2019) Study on the influence of soil microbial community 
on the long-term heavy metal pollution of different land use types and depth layers in 
mine. Ecotoxicology and Environmental Safety 170: 218-226. https://doi.org/10.1016/j. 
ecoeny.2018.11.136 

Zhuang P, McBride MB, Xia H, Li N, Li Z (2009) Health risk from heavy metals via consump- 
tion of food crops in the vicinity of Dabaoshan mine, South China. The Science of the 
Total Environment 407(5): 1551-1561. https://doi.org/10.1016/j.scitotenv.2008.10.061 


